0180  175  0EL8V  DISTRIBUTION  ANBLVSIS  OF  UINDOU  RANDOM-ACCESS  1/1 

ALGORITHMS  <U>  VIRGINIA  UNIV  CHARLOTTESVILLE  DEPT  OF 
ELECTRICAL  ENGINEERING  L  GEORG  I  AD  IS  ET  AL  29  APR  87 
UNCLASSIFIED  UVA/525415/EE87/105  N08014-87-K-0742  F/G  20/3  NL 


microcopy  resolution  test  chart 

NATIONAL  BORiAU  Of  STANDARDS  I96J  » 


AD- A 180  175 


► 


( 

I 


h 


\ 

\ 

I 


I 


OTIC  FILE  COK* 


A  Technical  Report 
Contract  No.  N00014-86-K-0742 
September  1,  1986  -  August  31,  1988 

DELAY  DISTRIBUTION  ANALYSIS  OF  WINDOW 
RANDOM-ACCESS  ALGORITHMS 

Submitted  to: 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  Virginia  22217-5000 

Attention:  Dr.  R.  N.  Madan 
Code  1114SE 

Submitted  by: 

L.  Georgiadis 

Research  Assistant  Professor 

M.  Paterakis 

Graduate  Research  Assistant 

P.  Kazakos 
Professor 


Report  No.  UVA/525415/EE87/105 
April  1987 

blSTRCBUTION  STATEMENT  A 
Approved  lor  public  releoael 
Distribution  Unlimited 


SCHOOL  OF  ENGINEERING  AND 
APPLIED  SCIENCE 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 


UNIVERSITY  OF  VIRGINIA 
CHARLOTTESVILLE,  VIRGINIA  22901 


A  Technical  Report 
Contract  No.  N00014-86-K-0742 
September  1,  1986  -  August  31,  1988 

DELAY  DISTRIBUTION  ANALYSIS  OF  WINDOW 
RANDOM-ACCESS  ALGORITHMS 

Submitted  to: 

Office  of  Naval  Research 
800  N .  Quincy  Street 
Arlington,  Virginia  22217-5000 

Attention:  Dr.  R.  N.  Madan 
Code  1114SE 

Submitted  by: 

L.  Georgiadis 

Research  Assistant  Professor 

M.  Paterakis 

Graduate  Research  Assistant 

P.  Kazakos 
Professor 


Department  of  Electrical  Engineering 


Report  No.  UVA/525415/EE87/105 


Copy  No. 


Unclassified 


SECuair 


la.  REP ORT  SECURITY  CLASSIFICATION 

UNCLASSI 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 

NONE  _ _ 


;  DISTRIBUTION/ AVAILABILITY  OF  KEFORT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION 
UNLIMITED 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

UVA/525415/EE87/105 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
University  of  Virginia 
Dept,  of  Electrical  Engr. 


6c.  ADDRESS  {City,  State,  and  ZIP  Code) 
Thornton  Hall 
Charlottesville,  VA  22901 


8a.  NAME  OF  FUNDING /SPONSORING 

organization  Dept,  of  Navy 
Office  of  Naval  Research 


8c.  ADDRESS  (City,  State,  and  ZIPCode) 

800  North  Quincy  Street 
Arlington,  VA  22217-5000 


6b.  OFFICE  SYMBOL 
(If  applicable) 


8b.  OFFICE  SYMBOL 
(If  applicable) 

N00014 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERlS) 


7a.  NAME  OF  MONITORING  ORGANIZATION 
Office  of  Naval  Research 
Resident  Representative _  _ 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

Joseph  Henry  Building,  Room  623 
2100  Pennsylvania  Avenue,  N.W. 
Washington,  DC  20037  _ 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

00014-87-K-0742 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


j  PROJECT 

TASK 

1  no 

NO. 

WORK  UNIT 
ACCESSION  NO 


11.  TiTLE  ( include  Security  Classification ) 

Delay  Distribution  Analysis  of  Window  Random-Access  Algorithms 


12  PERSONAL  AUTHOR(S) 

L.  Georgiadis,  M.  Paterakis,  P.  Kazakos _  _  _ 


13a.  TYPE  OF  REPORT  13b.  TIME  CCVEREO  14.  DATE  OF  REPORT  (fear,  Month.  Day)  i15.  PAGE  COUNT 

Technical _ from  9/1/86  TQ8/31/88  1  April  29,  1987  \  29 

16  SUPPLEMENTARY  notation 


COSATl  COOES 


GROUP  SU3-GROUP 


18.  SU8JECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


9  ABSTRACT  ( Continue  on  reverse  it  necessary  and  identify  by  block  number) 


A  method  for  delay  distribution  analysis  of  Window  Random-Access 
algorithms  is  presented.  The  window  size  is  allowed  to  vary  during  the 
operation  of  the  algorithm.  It  is  shown  that  the  quantities  of  interest 
in  the  computation  of  the  delay  distribution,  can  be  related  to  the 
solution  of  appropriate  infinite  systems  of  linear  equations.  Once  the 
constants  and  the  coefficients  of  the  unknowns  of  a  system  are 
determined,  bounds  on  the  solution  can  be  developed  by  applying 
previously  developed  methodologies.  The  method  is  applied  to  the  delay 
distribution  analysis  of  the  Capetanakis  Window  Random-Access  algorithm 
and  the  Fart-and-Try  algorithm,  both  under  binary  feedback. 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


20  US' 

•'  Jo’ 

ON  /  AVAilABUi  ’ 

'  •  A o S  *  a  A L  ' 

C3  ■  • 

c  AS 

t’ED'UNUMi't  0 

G  SAVF  AS  aPT 

□  O' C  St  PS 

22a  -A 

V  L  » j 

c  /NSiSLfc  % 

DIVIULAl 

i 

Dr. 

R. 

N.  Mad an 

1 

22b  IELLPhONE  {Include  Area  Code)  j  22c  OFHCc  Symbol 


(202)  696-4217 


OD  FORM  1473, 84  mar 


jj  Ai'J  eon.on  may  oe  uwo  until  «t«nausteo 
All  other  edition*  are  oosolete 


17  I  N00014 


SECURITY  CLASSIFICATION  Or  Tins  PAGE 

Unclassified 


1.  INTRODUCTION 


Window  Random-Access  algorithms  constitute  an  important  class  of  Multiple- 
Access  algorithms;  they  are  distributive  and  attain  high  throughputs  and  low  delays  by 
controlling  the  number  of  simultaneously  transmitting  users.  The  throughput  analysis  of 
algorithms  in  this  class,  is  a  relatively  easy  task.  The  delay  analysis,  however,  presents 
difficulties,  due  mainly  to  the  presence  of  variable  window  sizes  and  the  complicated 
state  space  that  some  of  these  algorithms  create.  This  restrictions  prohibits  the  applica¬ 
tion  of  results  from  standard  queueing  theory  in  the  delay  analysis. 

Many  attempts  for  the  delay  analysis  have  been  made.  In  [1],  the  class  of  algorithms 
with  constant  window  size  was  considered,  and  upper  bounds  on  the  expected  delays 
were  developed.  In  [2],  [3],  methods  for  the  computation  of  bounds  on  the  moments  of 
the  delays  were  presented.  A  method  for  the  computation  of  delay  distribution  for  con¬ 
stant  window  algorithms  appears  in  [4],  The  method  in  [4]  relies  on  a  clever  decomposi¬ 
tion  of  the  delay  process,  which  allows  the  application  of  results  from  standard  queueing 
theory.  It  is  not  possible,  however,  to  extend  the  techniques  in  [4]  to  the  case  of  variable 
window  size.  The  computation  of  the  delay  distribution  for  variable  window  size  algo¬ 
rithms  remains  an  open  problem.  One  possible  approach  is  to  compute  bounds  on  the 
moments  of  the  delays  as  in  [2]  or  [3],  which  can  then  be  used  for  an  approximate 
evaluation  of  the  delay  distribution.  This  approach,  however,  is  not  computationally  prac¬ 
tical. 

In  this  paper,  we  show  that  the  methodology  employed  in  [3],  can  be  extended  to 
provide  bounds  on  the  distribution  of  the  delays.  The  quantities  of  interest  are  related  to 
the  solution  of  a  denumerable  system  of  linear  equations.  Methods  for  the  computation 


of  the  constant  terms  and  the  coefficients  of  the  unknowns  of  the  system  are  developed. 
The  methodology  is  applied  to  the  delay  distribution  analysis  of  both  the  Capetanakis 
Window-Access  algorithm  with  binary  feedback  and  the  Part-and-try  algorithm  with 
binary  feedback.  It  can  also  be  applied  directly  to  other  Window  Random-Access  algo¬ 
rithms  with  different  feedback.  An  interesting  result  of  the  analysis  is  that  as  the  arrival 
rate  increases,  the  tales  of  the  distribution  become  longer,  but  the  median  grows  much 
slower  than  the  expected  delay. 

2.  MODEL  SPECIFICATION 

We  consider  a  single  slotted  channel  that  is  being  accessed  by  a  number  of  indepen¬ 
dent  packet  transmitting  users.  The  length  of  a  packet  is  equal  to  the  length  of  a  slot,  and 
packet  transmission  may  start  only  at  the  beginning  of  a  slot.  Simultaneous  transmission 
of  more  than  one  packets  in  the  same  slot,  results  in  complete  loss  of  the  information 
included  in  the  involved  packets.  The  latter  event  is  referred  to  as  a  "collision"  event.  At 
the  end  of  each  slot,  all  users  receive  a  feedback  that  provides  some  information  about 
the  chanell  activity  in  that  slot.  Common  types  of  feedback  are  the  binary  C-NC  (colli¬ 
sion  versus  noncollision)  feedback,  and  the  ternary  0-1 -C  (empty  versus  success  versus 
collision)  feedback.  To  resolve  the  collision,  the  users  follow  the  rules  of  a  Random- 
Access  algorithm.  The  algorithm  is  implemented  by  each  user  in  a  distributed  fashion, 
using  only  the  available  feedback.  The  cumulative  packet  generating  process  is  assumed 
to  be  Poisson  with  rate  X  packets  per  slot. 

We  assume  that  a  Window  Random-Access  algorithm  is  employed,  whose  basic 
operating  characteristics  are  the  following  (see  Figure  1):  Suppose  that  at  the  beginning 


of  slot  v  all  packet  that  arrived  before  time  ty<v  have  been  successfully  transmitted,  and 
there  is  no  information  concerning  the  packets  that  may  have  arrived  in  the  interval  [tv,v), 
(i.e.,  the  distribution  of  the  interarrival  times  of  the  packets  in  [ty,v)  is  the  same  as  the 
one  assumed  originally).  The  beginning  of  such  a  slot  v  is  called  a  Conflict  Resolution 
Point  (CRP).  The  time  difference  dy  =  v-ty  is  referred  to  as  the  lag  at  v.  In  slot  v,  the 
users  that  generated  packets  in  the  interval  [tv,  tv+tv),  where  tv  =  min(dv,A),  are  allowed 
to  transmit;  A  is  a  parameter  to  be  properly  chosen  for  throughput  maximization.  After  a 
random  number  of  slots  and  following  the  rules  of  the  algorithm,  another  CRP,  v ,  is 
reached,  with  a  corresponding  tv-  >  ty.  All  the  packets  that  have  been  generated  in  the 

interval  [ty  tv),  have  been  successfully  transmitted  in  the  interval  [v,v ).  The  intervals 

/ 

[v,v ),  [ty,  ty-Kty),  [tv,  tv)  are  called  conflict  resolution  interval,  transmitted  interval,  and 
resolved  interval,  respectively.  The  length  of  Tv,  is  called  the  window  size  at  time  v  . 
Clearly,  the  window  size  varies  with  time,  and  its  maximum  size  is  A.  Note  also,  that  the 
length  of  the  conflict  resolution  interval  is  one,  if  and  only  if  there  are  at  most  one  pack¬ 
ets  in  the  transmitted  interval. 

Algorithms  that  operate  as  described  above,  are  the  Capetanakis  Window  Random- 
Access  algorithm  [5]  and  the  Part-and-Try  algorithm,  [5],  [6],  under  either  binary  C-NC 
feedback,  or  ternary  0-1 -C  feedback. 

3.  STEADY  STATE  DELAY  DISTRIBUTION  ANALYSIS 

Let  packets  be  labeled  1,2,3,...  according  to  the  order  of  their  arrival  instants.  The 
delay  Dn  experienced  by  the  n-th  packet  is  defined  as  the  time  difference  between  its 
arrival  and  the  end  of  its  successful  transmission.  We  will  be  interested  in  evaluating  the 


steady  state  distribution  of  Da,  when  it  exists. 

Let  Vj ;  i>l  be  the  sequence  of  successive  CRPs  and  let  dj  be  the  lag  at  v;.  The 
sequence  dj ;  i>l  is  a  Markov  chain  with  state  space  F .  For  most  of  the  existing  Window 
Random- Access  algorithms,  F  is  a  denumerable  subset  of  the  interval  [1,°°).  Let  Tj=l, 
dj=l,  and  define  Ti+1,  as  the  first  slot  after  T;,  at  which  dT  =  1.  From  the  description  of 
the  algorithm  it  can  be  seen  that  the  induced  delay  process  probabilistically  restarts  itself 
at  the  beginning  of  each  slot  T;,  i=l,2,...  The  interval  [Tj,Ti+1)  will  be  referred  to  as  the  i- 
th  session.  Note  that  the  sessions  have  lengths  that  are  i.i.d.  random  variables. 

Let  Rj,  i=l,2,...,  denote  the  number  of  packets  successfully  transmitted  in  the  inter¬ 
val  (0,Tj];  (note  that  Rs  also  represents  the  number  of  packets  arrived  during  the  interval 
[0,  T-l),  since  T;  is  a  CRP  at  which  the  lag  is  unity).  Then,  C—R^j-Rj;  i>l,  is  the 
number  of  packets  successfully  transmitted  in  the  interval  (TjTj+1]  -  these  are  the  packets 
that  arrived  during  the  interval  [Tj-1,  Ti+1-1).  The  sequence  R;;  i>l,  is  a  renewal  process, 
since  C;,  i>l,  is  a  sequence  of  nonnegative  i.i.d.  random  variables.  Furthermore,  the 
delay  process  Dn,  n>l,  is  regenerative  with  respect  to  the  renewal  process  R;,  i>l,  with 
regeneration  cycle  C  x. 

Let 


1 

0 


ifo„a 

otherwise 


From  the  regenerative  theorem  [3],  we  conclude  that  if  C=E(Cj)  <  then 


P(x,r/x) 


:  The  probability  that  the  conflict  resolution  interval  has  length  x  and  the 
resolved  interval  has  length  r,  given  that  the  window  size  is  x. 


p(x/t)  :  The  probability  that  the  conflict  resolution  interval  has  length  x,  given 

that  the  window  length  is  x. 

hd  :  Number  of  slots  needed  to  reach  a  CRP  with  lag  1  given  that  the  current 

lag  is  equal  to  d,  d eF . 

kd(s)  :  Number  of  successfully  transmitted  packets  with  delay  less  than  s,  in  the 

interval  hd. 

d(s)  :  Number  of  successfully  transmitted  packets  with  delay  less  than  s  during 

a  conflict  resolution  interval,  given  that  the  window  size  is  x  and  the  lag  is 
d. 

nx(s)  :  Number  of  successfully  transmitted  packets  with  delay  less  than  s,  dur¬ 

ing  a  conflict  resolution  interval,  given  that  the  window  size  is  x  and  the 
lag  is  x.  That  is,  nT(s)=mr  T(s). 

Let  us  also  define. 


Note  that  by  definition, 


Also, 


Kd(s)  =  E(kd(s)) 
Mt>d(s)  =  ECm^Cs)) 
Nt(s)  =  E(nx(s)) 
Hd  =  E(hd) 


c, 

Ki(s)  =  E(£ln(s)) 

n=l 


7 


^-XH! 

Therefore,  the  determination  of  Kt(s)  and  Ht,  will  permit  the  computation  of  the  steady- 
state  distribution  of  the  delays. 

Consider  the  arrangement  of  Fig.  1.  The  delay  D  of  the  successfully  transmitted 
packet  1,  can  be  decomposed  as  follows: 


Therefore, 


D  —  0j+d— x+02 


D  <  s  iff  0,+0,  <  s-d+x 


But  0j+02  is  statistically  identical  to  the  delay  that  the  successfully  transmitted  packet 
experiences  if  the  transmitted  interval  is  x  and  the  lag  is  x.  The  last  observation  shows 
that  ntj  d(s)  is  identically  distributed  with  n^fs-d+x).  Observe  now  that 


kd(s)  = 


m^fs)  if  dv-  =  1 
ntj  d(s)+kd(s)  if  dv  =  d  *  1 


and  that 


dv  -  dy-S+Z ,  x  = 


dv  if  dv<A 
A  if  dv>A 


From  (4)  and  (5)  we  conclude  that 


Kd(s)  =  Mdd(s)  +  £Kd_r+x(s)p(x,r/d)  if  l<d<A,  deF 


Kd(s)  =  MAd(s)  +  £Kd_r+x(s)p(x,r/A)  if  d>A,  deF  (6 

r,x 

Since  m^s)  is  identicaly  distributed  with  nT(s-d+x),  equations  (6a)  and  (6b)  become. 


ii. 


»  1* 

I 


$ 

I 


§ 

$ 

I 


$ 

m 

3 

ft 


Kd(s)  =  Nd(s)  +  £Kd_r+x(s)p(x,r/d)  if  l<d<A,  deF 


(7a) 


(7b) 


Kd(s)  =  NA(s-d+A)  +  £kd_r+x(s)p(x,i7A)  if  d>A,  d£F 

r,x 

Equations  (7a),  (7b)  comprise  a  denumerable  system  of  linear  equati  'ns.  Of  interest  to  us 
is  the  element  K t(s)  of  a  particular  solution  of  this  system.  The  methodology  developed 
in  [3]  can  be  used  for  the  study  of  system  (7).  Note  that  the  coefficients  of  the  unknowns 
are  independent  of  s.  This  observation  represents  a  computational  advantage  when  the 
solution  to  (7)  is  approximated  by  the  solution  of  appropriate  finite  linear  system  of  equa¬ 
tions  [3],  In  this  case,  the  approximate  solution  can  be  represented  in  the  form, 

K(s)  =  (I-A)-1N(s) 

where  A  is  a  square  matrix  whose  elements  are  independent  of  s.  The  matrix  (I— A)-1  can 
be  computed  once,  and  then  used  for  the  computation  of  the  approximate  solution  for 
various  values  of  s. 

We  now  proceed  in  the  development  of  an  initial  upper  bound  on  the  solution  of 
system  (7).  Following  the  methodology  in  [3],  such  a  bound  will  be  the  sequence 
Kd(s)  =  Yu(s)d+^u(s),  if  Yu(s),  ^u(s),  can  be  determined  so  that  the  following  inequalities 


are  satisfied 


Kd°(s)  >  Nd(s)  +  £Kd°_r+x(s)p(x,r/d)  =  Kd(s)  if  l<d<A,  deF 


Kd°(s)  >  Nd(s-d+A)  +  XKd°.r+x(s)p(x,r/A)  -  Kd‘(s)  ifd>A.  deF 


Substituting  Kd(s)  in  the  right  hand  side  of  inequalities  (8),  it  can  be  easily  seen  that  if 


Kd‘(s)  =  Kd(s)  -  Nd(s)  +  Yu(s)(E(//d)  -  E(6/d)  -  (l+kd)e’>jl)  -  ;u(s)(l~Xd)e'^  if  l<d<A  (9a) 
Kd(s)  =  Kd°(s)  +  NA(s-d+A)  -  yu(s)(E(5/A)-E(//A))  if  d>A  (9b) 


Observe  now  that  NA(s)  is  an  increasing  function  of  s.  Therefore,  from  (9b)  we  conclude 
that, 


Kd'(s)  <  Kd°(s)  +  NA(s)  -  Yu(s)(E(5/A)-E(//A))  if  d>A  (10) 

From  (10)  we  conclude  that  if  E(//A)  <  E(5/A),  the  condition  for  stability  of  the  system, 


inequalities  (9b)  are  satisfied  if 


Yu(s)=' 


na(s> 

E(8/A)-E(//A) 


With  this  value  of  yu(s),  it  can  be  seen  that  inequalities  (9a)  are  satisfied  if 


where 


y(d)  = 


£u(s)  =  max{-yu(s),  sup  (\j/(d))} 

l<d<A 


Nd(s)  +  y„(s)[E(//d)  -  E(5/d)  -  (1+Xd)e"w] 


(1+Xd)e~ 


From  the  above  discussion  we  conclude  that  the  solution  to  system  (7)  satisfies  the  ine¬ 


qualities 


Kd(s)  <  Yu(s)d  +  Cu(s),  dzF 


where  yu(s),  Cu(s),  at®  given  by  equations  (11),  (12)  respectively.  The  uniqueness  of  the 
solution  is  guaranteed  by  the  same  techniques  as  in  [3].  If  we  use  a  similar  method  for  the 
development  of  a  lower  bound,  we  find  that 


where 


Y/(s)d  +  ^(s)  =  Kd(s),  deF 


Y/ (s)=0  and  ^(s)=  inf  {(Nd(s)/((1+Xd)e  )} 

l<d<A 


Bounds  on  Ht  are  given  in  [3],  formulas  (18),  (19). 

As  is  explained  in  [3],  the  bounds  (14)  can  be  used  to  further  improve  the  bounds  on 
K^s).  To  proceed  further,  however,  the  computation  of  the  quantities 
E(//d),  E(S/d),  p(x,r/d),  and  Nd(s)  is  necessary.  In  section  3,  we  present  a  method  for  the 
computation  of  Nd(s)  for  the  Capetanakis  dynamic  algorithm,  and  then  proceed  in  the 
computation  of  tight  upper  and  lower  bounds  on  the  quantities  of  interest.  In  section  4, 
we  present  a  method  for  the  computation  of  the  quantities  Nd(s)  and  p(x,r/d)  for  the 
Part-and-Try  algorithm,  and  then  we  develop  bounds  on  the  distribution  of  the  delays. 
Due  to  the  complicated  state  space  of  the  latter  algorithm,  the  development  of  tight 
bounds  for  high  input  rates  becomes  computationally  cumbersome. 

3.  THE  CAPETANAKIS  WINDOW  RANDOM-ACCESS  ALGORITHM  WITH  C-NC 
FEEDBACK 

A  complete  description  of  the  algorithmic  rules  can  be  found  in  [5].  In  this  algo¬ 
rithm,  the  resolved  interval  is  always  equal  to  the  window  size  i.e.,  x  =  8.  This  results  in 
a  significant  simplification  of  the  state  space  F ,  when  the  maximum  window  size  A  is  a 
rational  number.  The  restriction  of  A  to  the  rational  numbers  facilitates  the  development 
of  tight  bounds  for  the  delay  distribution  and  does  not  represent  any  disadvantage  in 
practice.  Moreover,  if  m  packets  are  involved  in  a  conflict,  the  conflict  resolution  process 
depends  only  on  m  and  not  on  the  window  length  or  the  generation  time  of  each  of  the 
packets.  The  last  property  facilitates  the  development  of  efficient  methods  for  the  com¬ 
putation  of  the  quantities  of  interest. 


Since  T  =  8,  we  have  that,  E(8/d)  =  d  and 


,  ...  f  p(x/d)  if  r=d 

P(u,d)={  0  otherwise 

Formulas  for  the  computation  of  E(//d),  p(x/d),  can  be  developed  by  following  the  rea¬ 
soning  in  [5].  In  the  next  section  we  develop  a  method  by  which  the  quantities  Nd(s)  can 
be  computed. 


3.1  The  Computation  of  Nd(s). 

Let  M  be  the  number  of  packets  in  the  window  d  (see  Fig  2).  Then, 

Nd(s)  =  E(nd(s))  =  E(E(nd(s)  /  M))  =  £  E(nd(s)  /  M  =  m)tU- -  (15) 

m=l  m! 

Let 


.1  if  the  delay  of  the  i-th  packet  in  d  is  at  most  s 
Ji-m(S)  =  i0  otherwise 


(16) 


The  "i-th  packet"  in  (16),  is  the  i-th  packet  in  a  random  enumeration  of  the  m  packets 
(randomly  chosen  packet).  Then 


m 

ni(s)  =  IU5)  (17) 

i=l 

Given  M,  the  generation  time  of  each  of  the  M  packets  is  uniformly  distributed  in  the 
window  d  and  independent  of  the  generation  time  of  the  rest  of  the  packets.  From  this 
observation  and  the  definition  in  (16)  we  conclude  that  the  random  variables 
Ji  m(s),  l<i<m.  are  identically  distributed  (although  not  independent).  Therefore,  from 


( 17)  we  conclude  that 


E(nd(s)  /  M  =  m)  =  mE(Jl  m(s)  /  M  =  m) 


(18) 


Let  0j,  (02),  be  the  delay  of  the  first  randomly  chosen  packet  before  (after)  the  initializa¬ 


tion  of  the  collision  resolution  process  i.e.  until  (after)  time  v  (Fig.  2).  Since  0j  is  uni¬ 


formly  distributed  in  [0,  d),  we  have  that 


Ed!  m(s)  /  M  =  m)=  j  Ed!  m(s)  /  M  =  m^  =  0  )  d0 
d  o 


Since  the  conflict  resolution  process  is  independent  of  the  packet  generation  time  in  a 


window,  we  conclude  that  given  m,  the  random  variables  0j,  02,  are  independent.  There- 


Ed^s)  /  M  =  m,  0j  =  0)  =  p(0j+02  <  s  /  M  =  m,©!  =  0 )  =  p(02  <  s-0  /  M  =  m)  (20) 


E(Ji  m(s)  /  M  =  m)  =  —  J  p(02  <  s-0  /  M  =  m)  d0 
d  o 


Observe  now,  that  02  takes  only  positive  integer  values.  Let 


Pqm)  =  P(02  =  q  /  M  «  m) ;  q=l,2,3,... 


then, 1,2 


p(02  2  s-0  /  M  =  m)  =  £ 

q=l 

Substituting  (22)  in  (21),  we  conclude  that 


1  ls-ei 

E(Jlm(s)/M=m)  =  -/  IPq(m)d0 

d  0  q=l 


[a  |  denotes  the  integer  pan  of  a. 


1  We  adopt  the  notation  £a.,=0  if  n>m 


ft 


WWW.. 
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Let  e=d-s+lsj-  Then,  since  P ^  is  independent  of  9,  (23)  can  be  written  as  follows: 

I  (s|  [s|-l  |5)-le|  -  lsj-lej-l 

E(Jlm(S)/M=m)  =  -((s-tsl)5;P(!m)+  IPq(m)+-.+  I  P‘m)  +  0HeJ)  I  P'm))  (24) 

q=l  q=l  q=l  q=l 

From  (15),  (18)  and  (24),  we  observe  that  for  the  computation  of  Nd(s),  only  the  quanti¬ 
ties  Pqm)  are  needed.  In  Appendix  1,  we  provide  recursive  formulas  for  the  computation 
of  the  quantities  Pqffl).  It  is  worth  noticing  that  formulas  (15),  (18)  and  (24),  are  valid  for 
any  window  Random- Access  algorithm  that  has  the  properties  described  in  the  first  para¬ 
graph  of  section  3. 


3.-2  Development  of  bounds  on  Kt(s),  H, 

For  the  development  of  bounds  on  the  delays  we  chose  A=2.5  As  a  result,  the  state 
space  F  becomes  simply, 

F={1,  1.5,2.,  2.5,....} 

For  the  algorithm  considered  in  this  section,  the  maximum  throughput  is  achieved  for 
A=2.67  [5].  The  reduction  in  throughput  due  to  the  choice  A=2.5  is  insignificant  (less 
than  .1%). 

Since  8=rc,  equations  (7a),  (7b),  become: 


Kd(s)  =  Nd(s)  +  XK*(s)p(x/d)  if  <j=l,  1.5,  2,  2.5 

X 

%m\ 

Kd(s)  =  NA(s-d+A)  +  £Kd_A+)l(s)p(x/A)  if  d=3,  3.5,  4,... 


(25a) 

(25b) 


X 

For  the  development  of  bounds  on  K,(s)  we  followed  the  method  of  truncation  of  the 


infinite  system  (25)  [3],  Specifically,  in  system  (25),  we  replaced  the  unknowns 
Kd(s),  for  d>40  with  the  upper  bounds  in  (14a).  The  substitution  results  in  a  finite  system 
of  equations  whose  solution  is  an  upper  bound  to  the  solution  of  (25)for 


d=l,  1.5,  2,. ..,40..  For  the  development  of  lower  bounds,  we  replaced  the  unknowns 
Kd(s),  d>40  with  the  lower  bound  in  (14b).  The  same  methodology  is  employed  for  the 
development  of  bounds  on  Hr  The  resulting  upper  and  lower  bounds  on  the  distribution 
of  the  delays  differ  by  .01  in  the  worst  case.  In  figure  3,  we  provide  the  distribution  of 
the  delays  for  various  values  of  the  arrival  rate.  An  important  observation  is  that  as  the 
arrival  rate  increases,  the  tales  of  the  distribution  become  longer,  but  the  median  grows 
much  slower  than  the  expected  delays. 

4.  THE  PART-AND-TRY  ALGORITHM  WITH  C-NC  FEEDBACK 

A  detailed  description  of  the  algorithm  can  be  found  in  [6].  The  algorithm  has 
throughput  .45,  operates  under  binary  feedback  and  it  is  not  blocked  in  the  presence  of 
feedback  errors.  The  techniques  used  in  the  present  section,  can  be  easily  applied  to  the 
analysis  of  the  algorithm  under  ternary  feedback. 


4. 1  The  computation  of  >L 


Let  M  be  the  number  of  packets  in  the  window  d  (see  Fig  4.).  As  in  section  3.1,  for¬ 
mula  (15)  holds.  Since  not  all  the  packets  in  a  window  are  successfully  transmitted,  how¬ 
ever,  we  need  to  modify  the  definition  of  Jj  m(s)  in  (16).  Let 

{1  if  the  i-th  randomly  chosen  packet  in  d  is  sucessfully 

transmitted  and  its  delay  is  at  most  s.  (26) 

0  otherwise 


then,  as  in  (18), 


and 


E(nd(s)  /  M=m)  =  mE(7  lm(s)  /  M=m) 


(27) 


1 

E(7  j  m(s)  /  M=m)  =  -  J  E(7 1  m(s)  /  M=m,  91=6)d6  (28) 

d  o 

Again,  02  takes  integer  values.  However,  0t  and  02  are  not  independent.  To  proceed 
further,  we  define  the  events 

A  i>m(s)={  the  i-th  randomly  chosen  packet  is  successfully  transmitted  and  02<s.} 
aim(q)  =  (the  i-th  randomly  chosen  packet  is  successfully  transmitted  and  92=q,  q=l,2,...} 

Since  0  takes  only  integer  values,  we  have  that 

UJ 

'4i,m(s)=  (29) 

q=l 

From  (28)  and  (29)  we  conclude  that 

d 

E(7 !  m(s)  /  M=m)  =  -J  p (A  ,  m(s-9)  /  M=m,  O^OMO  (30) 

d  o 

l  d  ls_^l 

=  t j  Z,  p(<*ijm(q) ' M=m’  ei=0)de 

d0  q=l 

Observe  now  that  p (a  l  m(q)  /  M=m,  9j=9)  depends  on  0  through  the  ratio  <>=9/d.  Let  us 
define 

0 

fi(q,m,<J>)  =  p(a s  m(q)  /  M=m,  O^c^d) ;  <>=—  e  [0,1)  (31) 

d 

then,  the  following  equations  hold: 


f^O.m.i!))  =  0,  for  m=l,2,...,  and  <>£[0,1 ) 
fL(  1  ,m,<»  =  0,  for  m=2,3,..,  and  <>£[0,1 ) 

ot'hfxL ; 


(32a) 

(32b) 

(32c) 


lit— I  , 

Sf1(q-lj+l,2<t>-l)  m“  2Hm_1)) 
j=o  -  J  - 


.5  <<(><  1 


fi(q,m,<!>)  = 

"*2  lf1(q-2,m,24.)2-<m-1)  +  f1(q-2,m-l,2<i))(m-l)2(-<m-1))  Q^<5 


Formulas  (32a),  (32b),  (32c),  are  obvious.  Let  us  explain  formula  (32d).  Assume  that 
.5<4»<1.  Then,  the  packet  under  consideration  (packet  1  in  fig  4),  lies  in  the  left  hand  half 
(l.h.h.)  of  the  window.  Let  02  be  the  delay  of  packet  1,  after  thcfirst  collision.  Then, 


02  =  02+l 


The  location  of  packet  1  in  the  l.h.h.  of  the  window,  is, 


0j  =  0r-d/2 


Therefore,  the  new  ratio  becomes 


<b  = - =  20-1 

d/2 


Let  Gj  m_1  be  the  set 


G  j  m_j={j  of  the  m-1  packets  (other  than  packet  1),  are  located  in  the  l.h.h.  of  the  window. } 


Note  that  p (G 


-  f 


-1  --(m-1) 


2  1  \  From  (33),  (34)  and  (35),  we  derive  formula  (32d) 


for  .5<4>^1  by  conditioning  on  Gj  If  0<<}><.5,  packet  1  is  located  in  the  right  hand  half 
(r.h.h.)  of  the  window,  In  this  case,  we  note  that  if  in  the  l.h.h.  of  the  window  there  are 
more  than  one  packets,  then  packet  1  is  not  transmitted  during  the  collision  resolution 
process.  By  conditioning  again  on  the  Gj  and  taking  into  account  the  last  observa¬ 
tion,  we  derive  equation  (32d)  for  0<<J><.5. 


■|Vtl *>•>  !■. V  '  VA  y  .  IV  Vq 


v.v 


From  equations  (32)  we  conclude  by  induction  the  following  property  for  the  func¬ 


tion  fjCq.m.^y. 

Property  1.  For  fixed  q  and  m,  the  function  fjCq.m.O)  is  simple  (i.e.,  it  takes  a  finite 
number  of  values),  and  left  continuous.  The  jumps  of  the  function  occur  at  the  points 
0,2-(q-1), ....  k2“(q_l), ....  l-2~(q_1). 

Taking  into  account  Property  1,  we  can  compute  E(/j  m(s)),  using  formula  (30). 

4.2  The  Computation  of  p(x,r/d). 

The  method  of  conditioning  on  the  number  of  packets  in  a  window  (applied  in  sec¬ 
tions  3.1  and  4.1),  does  not  seem  to  lead  to  easily  computable  recursive  formulas.  In  this 
section,  we  present  an  alternative  methodology  that  results  in  simple  recursive  formulas 
for  the  computation  of  p(x,r/d),  and  provides  insight  into  the  structure  of  these  probabili¬ 
ties. 

Observe  first  that  p(x,r/d)  depends  on  r  through  the  ratio  s=r/d.  Let 

f2(x,s,d)  =  p(/  =x,  8=sd,  M>2  /  d) 

Conditioning  on  the  events  {M=0},  {M=l}  and  {M>2},  and  observing  that  if  M=0,1  then 
/=1  and  8=d,  we  conclude  that 

p(x,r/d)  =  8K(x-l)8K((s/d)-l)e_(Xd)(l+tai)  +  f2(x,s,d)  (36) 

where 

5K(y) 

The  function  f-,(x,s,d)  satisfies  the  following  recursive  formulas 


fl  if  y=0 
jo  otherwise 


f2(l,s,d)  =  f2(2,s,d)  =  0,  0<s<l 


(37a) 


f2(x-l,2s,d/2)  0<s<5 

f2(x-2,2s-l,d/2)e~(X^2)(l+Xd/2)  -5<s<l 
+  8K(x-3)SK(s-l)e~(Xd)(Xd/2)2 


(37b) 


We  derive  equation  (37a)  for  0<s<  5  The  case  .5<s<l  can  be  derived  by  a  similar  reason- 

*  / 

ing.  Let  M  be  the  number  of  packets  in  the  l.h.h  of  the  window.  Let  /  be  the  number  of 

slots  during  a  collision  resolution  process  after  the  first  collision.  Then, 

p(i  =x,  5=sd,  M>2  /  d)  =  p(/'+l=x,  §=sd,  M>2,  M>2  /  d)  +  p (/=x,  S=sd,  M>2,  MSI  /  d) 

But  if  M<1,  then  s>.5.  Therefore  p(/=x,  8=sd,  M>2,  M<1  /  d)  =  0  if  0.<s<5.  It  remains 
to  observe  that  p(/  +l=x,  8=sd,  M>2,  M  >2  /  d)  =  p (/  =x-l,  8=2(sd/2),  M  >2  /  d),  and  that 
if  M  >2,  the  collision  resolution  process  after  the  first  collision  is  statistically  identical  to 
a  collision  resolution  process  that  started  with  a  window  of  size  d/2. 

From  (37a),  (37b)  we  conclude  by  induction  that  the  function  f2(x,s,d)  has  the  fol¬ 
lowing  property: 

Property  2  For  fixed  x>3,  and  for  any  d,  s  takes  a  finite  number  rx  of  values.  The 
sequence  rx  satisfies  the  following  recursions:  r3=r4=l,  rx=rx_i+r*_2->  x>5.  Let  A  x  be  the 
set  of  values  of  s  for  given  x.  Then,3 

A3={1},  A4={  1/2},  Ax  =  (.5Ax_,)u(.5>4x_2+.5);  x>5 

The  values  of  f2(x,s,d)  for  various  x  and  s,  can  be  computed  from  (37b).  The 


1  We  use  (he  nouuon  +b  =  (y  :  y=ax+b,  xe4  } 


probabilities  p(x,r/d),  are  easily  computed  from  (36). 

4.3  The  Development  of  bounds  on  K^sl.H,. 

The  state  space  F  for  the  Part-and-Try  algorithm  is  a  dense  subset  of  [1,  °°).  This 
property  complicates  the  development  of  finite  systems  of  linear  equations  whose  solu¬ 
tion  will  provide  upper  or  lower  bounds  to  the  quantities  K,(s),  H{.  Taking  advantage  of 
the  structure  of  the  probabilities  p(x,r/d)  (see  section  4.2),  however,  we  can  proceed  as 
follow  s:  Let  B  be  a  finite  subset  of  F ,  that  includes  the  state  d=l.  We  develop  the  follow¬ 
ing  finite  system  of  linear  equations: 

Yd  =  Nd(s)  +  ]£  (Yu(s)(d-r+x)  +  ^u(s))p(x,r/d)  +  £  Yd_f+lp(x,r/d)  if  l<d<A,  de8 

(d-r+*)tfl-  (d-r+*)e8  (38a) 

Wl 

Yd  =  NA(s-d+A)  +  £  (Yu(s)(d-r+x)  +  £u(s))p(x,r/d)  +  £  Yd_f+Ip(x,r/A)  if  d>A.  de£ 

(d-r+x)efl*  (d-r +x)cB 

yu(s)  and  £u(s)  are  determined  from  (11),  (12).  Due  to  property  2  of  section  4.2,  it  is  sim¬ 
ple  to  find  for  a  given  defi ,  the  values  of  x  and  r  such  that  (d-r+x)  e  B .  The  summation 
over  the  infinite  set  Bc,  can  be  computed  in  terms  of  E(//d),  E(o/d),  yu(s),  £u(s),  and  the 
probabilities  p(x,r/d);  (d-r+x)  e  B .  The  solution  Yd;  de5  of  system  (38),  is  an  upper 
bound  to  the  solution  Kd(s);  deB  of  system  (7)  [3].  Since  leB ,  we  can  determine  a  bound 
on  Kj(s).  Similarly,  lower  bounds  on  K^s),  and  upper  and  lower  bounds  on  Hj  can  be 

developed.  For  the  computations  we  used  the  set  B  ={l,  1.125,...,  l+(k/8) . 9}. The 

resulting  bounds  on  the  delay  distribution  are  presented  in  Fig  5  for  X  =  .1,  .2,  .3.  For 
higher  arrival  rates,  the  bounds  are  not  tight  (see  Table  1)  and  although  they  can  be 
improved  by  enlarging  the  set  B,  the  computations  become  cumbersome. 


1 


5.  CONCLUSIONS 

We  developed  a  method  for  the  computation  of  bounds  on  the  delay  distribution  of 
Window  Random- Access  algorithms.  The  method  has  been  applied  to  the  delay  distribu¬ 
tion  analysis  of  the  Capetanakis  Window  Random-Access  algorithm  and  the  Part-and- 
Try  algorithm  both  under  binary  C-NC  feedback.  The  bounds  developed  for  the 
Capetanakis  algorithm,  are  tight  for  all  arrival  rates  within  the  stability  region  of  the 
algorithm.  For  the  Part-and-Try  algorithm,  however,  the  bounds  are  satisfactory  for  rela¬ 
tively  low  arrival  rates.  The  computational  difficulty  in  obtaining  tight  bounds  for  the 
latter  algorithm,  is  due  to  its  complicated  state  space.  The  techniques  can  be  easily 
applied  to  other  Window  Random-Access  algorithms  whose  operating  characteristics  are 
as  described  in  Section  2. 
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APPENDIX 


Here,  we  provide  recursive  formulas  for  the  computation  of  P^ .  Clearly, 


p(m)  =  0-q<0  P^^i1  1 

q  ’ q  ’  q  [0  otherwi: 


q  ^  [^0  otherwise 

Let  /n  be  the  number  of  slots  needed  for  the  resolution  of  multiplicity  n>0  conflict.  Then, 


from  the  operation  of  the  algorithm  we  conclude  that 


1  f  m-l]  1 


r  1  f  m-l]  1 

p(n+l)  with  prob.  —  - 

rq-l  2  L  n  J  'J®-1 

p(m)  _ .  *• 

1 2  With  prob.  -rf1"”1]  p(/n  =  0 


(A.2) 


The  upper  part  of  (A.2)  is  derived  by  considering  the  event  that  the  packet  under  con¬ 
sideration,  together  with  n  of  the  rest  m-l  packets  retransmit  immediately  after  the  initial 
collision.  The  lower  part  of  (A.2)  is  derived  by  considering  the  event  that  the  packet 
under  consideration  does  not  transmit  immediately,  while  the  n  of  the  rest  m-l  packets 
retransmit  immediately  after  the  initial  collision  and  it  takes  /  number  of  slots  to  resolve 
a  collision  of  multiplicity  n.  The  probabilities  p(/n=/ )  can  be  computed  by  similar  rea¬ 
soning.  Averaging  in  (A.2),  we  finally  have  the  following  recursive  formulas  for  P*m)  for 


in— 1  r  / 

p(m)  _  ~-m  ^-i  p(n+l)  rn— 1  . 

* q  ^  l  n  2 


m-l  q-2  r 

-m  r,  y  ml  ^(m-n) 

■  L  L[  n  \  *q-W 

D=0  /=1 


P('„=') 
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Table  1 


Bounds  on  the  delay  distribution  for  the  Part-and-Try  algorithm 
with  binary  C-NC  feedback.,  for  A.  =  0.4. 
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UNIVERSITY  OF  VIRGINIA 
School  of  Engineering  and  Applied  Science 

The  University  of  Virginia’s  School  of  Engineering  and  Applied  Science  has  an  undergraduate 
enrollment  of  approximately  1,500  students  with  a  graduate  enrollment  of  approximately  560.  There 
are  150  faculty  members,  a  majority  of  whom  conduct  research  in  addition  to  teaching. 

Research  is  a  vital  part  of  the  educational  program  and  interests  parallel  academic  specialties. 
These  range  from  the  classical  engineering  disciplines  of  Chemical,  Civil,  Electrical,  and  Mechanical 
and  Aerospace  to  newer,  more  specialized  fields  of  Biomedical  Engineering,  Systems  Engineering, 
Materials  Science.  Nuclear  Engineering  and  Engineering  Physics,  Applied  Mathematics  and  Computer 
Science.  Within  these  disciplines  there  are  well  equipped  laboratories  for  conducting  highly  specialized 
research.  All  departments  offer  the  doctorate;  Biomedical  and  Materials  Science  grant  only  graduate 
degrees.  In  addition,  courses  in  the  humanities  are  offered  within  the  School. 

The  University  of  Virginia  (which  includes  approximately  2,000  faculty  and  a  total  of  full-time 
student  enrollment  of  about  1 6,400),  also  offers  professional  degrees  under  the  schools  of  Architecture, 
Law,  Medicine,  Nursing,  Commerce,  Business  Administration,  and  Education.  In  addition,  the  College 
of  Arts  and  Sciences  houses  departments  of  Mathematics,  Physics,  Chemistry  and  others  relevant 
to  the  engineering  research  program.  The  School  of  Engineering  and  Applied  Science  is  an  integral 
part  of  this  University  community  which  provides  opportunities  for  interdisciplinary  work  in  pursuit 
of  the  basic  goals  of  education,  research,  and  public  service. 


